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ABSTRACT 

> 

' Context. We present high resolution VLT UVES and low resolution FORS optical spectroscopy of supernova 2006aj and its host 

£f>^ ' galaxy, associated with the nearby (z = 0.03342) gamma-ray burst GRB 060218. This host galaxy is a unique case, as it is one of the 

few nearby GRB host galaxies known, and it is only the second time high resolution spectra have been taken of a nearby GRB host 
^ . galaxy (after GRB 980425). 

Aims. The resolution, wavelength range and S/N of the UVES spectrum combined with low resolution FORS spectra allow a detailed 
£^ ' analysis of the circumburst and host galaxy environments. 

Methods. We analyse the emission and absorption lines in the spectrum, combining the high resolution UVES spectrum with low 
resolution FORS spectra and find the metallicity and chemical abundances in the host. We probe the geometry of the host by studying 
the emission line profiles. 

Oh' Results. Our spectral analysis shows that the star forming region in the host is metal poor with 12 + log(0/H) = 7.54^ j 7 , (~ 0.07 Zq), 

placing it among the most metal deficient subset of emission-line galaxies. It is also the lowest metallicity found so far for a GRB 
host from an emission-line analysis. Given the stellar mass of the galaxy of ~ 10 7 Mq and the SFRho- = 0.065 ± 0.005 IviQyr -1 , 
the high specific star formation rate indicates an age for the galaxy of less than ~ 200 Myr. The brightest emission lines are clearly 
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asymmetric and are well fit by two Gaussian components separated by ~ 22 km s -1 . We detect two discrete Nal and Call absorption 

components at the same redshifts as the emission components. We tentatively interpret the two components as arising from two 
different starforming regions in the host, but high resolution imaging is necessary to confirm this. 

Key words. Gamma rays: bursts - galaxies: high redshift, abundances - cosmology: observations 

1. Introduction Universe (see e.g. Jakobsson et al. 2005). The significant in- 
crease in detection rate and localisation of long GRBs through 

Long-duration gamma-ray bursts (GRBs) are widely accepted ^ successful opera tion of Swift in principle permits the study 

to be related to core-collapse supernovae: clear supernova sig- of a j and uniforml selected le of GRB host laxies 

natures are seen in the afterglow spectra of low redshift GRBs Howeve] . t0 date the le of spectroscopica u y studied GRB 

£.g Stanek et al. 2003, Hjorth et al. 2003; Pian et al. 2006). hofjt laxies h smaU ( „ 3Q cafjes) and suffer from several 

Dedicated surveys of GRB hosts suggest that GRBs occur pref- selection biase ^ dug tQ theif faintness 

erentially in low mass, subluminous, blue star-forming galax- _ . , , . , . 

i Viu » i ™o -c u. * i man t Si >u » One of the key properties needed to understand GRB pro- 

les (e.g. Chary et al. 2002; Fruchter et al. 1999; Le Floe h et . , , . 3 F " . ........ . 

i inmi rr U V>r>D a. i » a uvu ■ u, e gemtors and their environments is the metallicity distribution 

al. 2003). The GRBs are often located within UV-bnght parts of f „ XT „ * „~„~ 

a. ■ u . u » e , , , /t)1 ° J oaat (e.g. Langer & Norman 2006; Yoon, Langer & Norman 2006): 

their hosts, where star formation takes place (Bloom et al. 2002; \ 6 , „ , , , ^^t? • , , 

c u, » i onn/;i j u \ u * * a the popular collapsar model for long GRBs requires low metal- 

Fruchter et al. 2006), and are shown to be more concentrated ,. . F F . F _ , ..... „ ™„ , , . 

■ , • , . . • e ■ , . ■ , licity progenitors. Metalhcities for GRB hosts can be determined 

towards the brightest regions of their hosts than are, in general, , , r , , , , , 

rc t i onn/;\ through afterglow spectroscopy and through host galaxy spec- 
core-collapse supernovae (Fruchter et al. 2006). & ., & ... ^TTT , 1 T 1 

r^nn u » 1 • » 1 * a lL u ,u • 1 • troscopy. Absorption lines of HI and heavy metals have pro- 

GRB host galaxies are not selected through their luminosities ., ,,. . " , „„„ . , ,. . , , , .„ . , 

. ,f 1 u ..I. £ 4. ..I. 4. i^t>t> i. u a ^ ^ a vided metalhcities along GRB sight-lines in the redshift interval 

or colours, but merely by the fact that a GRB has been detected. & b 

This could potentially provide an unbiased sample of starform- z 

ing galaxies, which may be used to study star formation in the Host g alax Y spectroscopy can provide metalhcities for the 

galaxy as a whole, but relies on the detection of the nebular 

Send offprint requests to: kwrsema@science.uva.nl emission lines, which is difficult at z > 1. There are only a few 

* Based on observations obtained at the ESO VLT under ESO pro- GRB host galaxies that are bright enough to permit a direct abun- 

gramme 076.A-0737(A). dance analysis through an electron temperature (T e ) determi- 
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nation (GRBs 980425, 020903, 031203; Prochaska et al. 2004, 
Hammer et al. 2006). It is only possible to study the metallicities 
of higher redshift and fainter hosts through secondary metallicity 
indicators using bright (nebular) lines (e.g. the R23 method), us- 
ing empirical correlations between the fluxes of certain emission 
lines and metallicity. 

On the 18th of February 2006 a bright, nearby GRB was dis- 
covered by Swift. The proximity (z ~ 0.0334, Mirabal & Halpern 
2006, the second closest GRB) of this GRB triggered a large 
follow-up campaign by several groups, which provided a unique 
opportunity to unveil the nature of a faint nearby galaxy associ- 
ated with a GRB event. GRB 060218 was found to be unusually 
long in duration (Tgo = 2100 ± 100 s, Campana et al. 2006) 
and of relatively low luminosity (Ey^o - 6.2 + 0.3 x 10 49 erg, 
Soderberg et al. 2006). Its prompt emission was soft, placing it 
in the class of X-ray flashes (XRFs). A bright supernova (des- 
ignated SN 2006aj) was clearly associated with this event which 
was studied with very high spectral and time resolution over a 
wide range of wavelengths, from X-rays to radio (e.g. Campana 
et al. 2006; Pian et al. 2006; Mazzali et al. 2006; Soderberg et 
al. 2006). The host was found to be a small galaxy with an irreg- 
ular morphology. 

In this paper we study the host of GRB 060218 through 
a high resolution VLT Ultraviolet and Visual Echelle 
Spectrograph (UVES; Dekker et al. 2000) spectrum, taken 
around the peak magnitude of the supernova. This spectrum 
shows a variety of well resolved emission lines, associated with 
ionized gas in the starforming regions in the host. The neutral 
gas of the ISM is probed by the detection of a few narrow ab- 
sorption lines. Low resolution FORS spectra (described in Pian 
et al. 2006) are used to study the fluxes of emission lines that fall 
outside the spectral coverage of the UVES spectrum. 

Several papers have already been published on GRB 060218, 
its associated supernova, and the host. However, on the metallic- 
ity of the host there is a wide range of reported values. Modjaz et 
al. (2006) report a metallicity of 0.15 Z , Mirabal et al. (2006) 
derive 0.46 Zq, while Sollerman et al. (2006) mention that the 
abundance is below Solar, but that the exact value is uncon- 
strained from the strong emission lines. These values are derived 
from different spectra, but the spread is largely due to internal 
scatter in empirical, secondary calibrators that are used, and the 
limited range of metallicities and diagnostic emission line ratios 
for which these secondary calibrators are valid (see e.g. Ellison 
& Kewley 2006). Since the derived metallicities are frequently 
used to draw conclusions on important issues such as progenitors 
(e.g. Sollerman et al. 2005) or the rate of nearby GRBs (Stanek 
et al. 2006; Wolf & Podsiadlowski 2006), it is clearly of great 
interest to pin down these uncertainties. 

The paper is organized as follows: In §2 we describe the ob- 
servations. In §3 we calculate the metallicity of the dominant 
starforming region(s) in the host galaxy and relative abundances 
of Ne and N, as well as the star formation rate from optical and 
radio fluxes. In §4 we analyse the discrete velocity components 
in the host through emission and absorption lines. We discuss the 
use of metallicity values found from secondary metallicity cali- 
brators in §5. In §6 we discuss the implications that a future de- 
tection of Wolf-Rayet star signatures in the host of GRB 060218 
may have on single star GRB progenitor models. 

Throughout this paper we use the cosmological parameters 
H = 70 km s '/Mpc, Q. M = 0.3 and Q A = 0.7. 



2. Observations 

GRB 0602 18/SN2006aj was observed with ESO VLT Kueyen 
(UT 2) on March 4, 2006, roughly at the time of maximum light 
of the supernova (around March 1, e.g. Sollerman et al. 2006). 
The UVES observation started at 00:30 UT, for a total exposure 
time of 2100 seconds. The magnitude of SN2006aj was V ~ 
17.6 at the time of observation. The airmass was high, averaging 
2.5, and the seeing measured from the 2D spectrum is 1 . 1 arcsec. 
The UVES-setup spectral resolution is R ~ 46000 (FWHM 
6.5 km s _1 ). The slit width was set at 1 arcsec, corresponding 
to about 7 kpc physical size at the distance of the GRB. The 
exposure was performed at the parallactic angle. The spectrum 
(wavelength ranges 3285 - 4527 A, 4621 - 5598 A and 5676 - 
6651 A) was reduced in the standard fashion using MIDAS and 
IRAF routines. 

The L.A. Cosmic program (specifically the lacosspec rou- 
tine, Van Dokkum 2001) was used to remove both point- and 
irregular shaped cosmic ray hits from the 2D spectrum before 
extraction. The spectrum was dispersion corrected, and flux cal- 
ibrated by archived response functions. An air to vacuum con- 
version and a heliocentric correction (+28 km s~') were applied 
to the spectrum. The resulting spectrum was compared with a 
(quasi)-simultaneous FORS spectrum taken 20 minutes after the 
UVES spectrum, which is flux calibrated through a standard star 
observation and through simultaneous photometry in B, V, R at 
VLT FORS2 (Pian et al. 2006), using the magnitude to flux con- 
versions from Fukugita, Shimasaku & Ichikawa (1995). To the 
UVES spectrum a multiplication factor of 1.625 was applied to 
match the well calibrated FORS flux values (compensating for 
the slit loss between the FORS and UVES spectra). 

We find a good match between the UVES and FORS spectra 
and photometry in the red end (above A ~4500 A), as shown in 
Figure Q] In the blue end (A < 4500 A) the UVES continuum 
flux is slightly higher than the FORS continuum flux. The prime 
reason is likely the very high airmass at which both spectra have 
been taken, making the flux calibration of both the UVES and 
FORS spectrum at the blue end more uncertain. We decide to 
not alter the flux calibration, but warn that the fluxes of emission 
lines below ~4500 A have a small additional uncertainty. This 
uncertainty does not significantly affect the metallicity results of 
this paper, because the electron temperature uncertainty is dom- 
inated by the uncertainty in the flux of [OIII] /14364 from the 
FORS spectra. The host galaxy emission line flux ratios found 
in the UVES spectrum agree, within the errors, with those found 
from the combination of the FORS spectra. 

A Galactic extinction correction was performed using E(B - 
V)mw = 0.142 (Schlegel, Finkbeiner & Davis 1998), assuming a 
Galactic extinction law A^/Ay expressed as Ry = Ay/E(B - V) 
(Cardelli et al. 1989), and Ry — 3.1. This value is slightly higher 
than the extinction derived by Guenther et al. (2006), who used 
the Galactic sodium lines to find E(B-V)mw = 0.127. Given that 
the systematic error in the conversion of the equivalent widths of 
Na to E(B - V) is poorly known, we choose to use E(B - V)mw 
= 0.142 mag. 

A selection of lines from the spectrum and the associated 
error spectra are shown in Figures [2] and [5] Although the spec- 
trum is dominated by the SN spectrum (V ~ 17.6 vs the host 
V ~ 20.2), several bright emission lines from the host galaxy are 
detected and resolved in the UVES spectrum (Fig. [2] Table [3). 
On the other hand, the detected absorption lines, from both the 
host galaxy and Milky Way, are narrower than the UVES resolv- 
ing power. The highest signal-to-noise emission lines provide 
a heliocentric mean redshift z = 0.03342(2) (Pian et al. 2006), 
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Fig. 1. The UVES spectrum shown in black, with overlaid in red 
the low resolution FORS spectrum. For presentation purposes 
the UVES spectrum has been smoothed to the pixelscale of the 
FORS spectrum. The points denote B, V, R VLT FORS photom- 
etry from the same night (Pian et al. 2006). The widths of the 
broadband filters are not shown. 



where we adopt a conservative error on the redshift due to poorly 
known systematic effects of the spectrograph (e.g. the centering 
of the object in the slit). 

The analysis of the emission lines was done using the IRAF 
software packages, mainly using the splot routines. As the emis- 
sion lines are clearly asymmetric (see §4.1), fluxes are measured 
using the numerical integration method (<? in splot), and not the 
standard Gaussian fitting. The errors in the line fluxes are gener- 
ally dominated by the uncertainty in the continuum level and are 
given at the l<x level. 



3. General host galaxy properties 

3.1. Metallicity 
3.1.1. Reddening 

Before we can use the emission lines in the FORS and UVES 
spectra to derive the host properties, we need to correct for the 
dust extinction intrinsic to the host. We use the Balmer line 
fluxes to measure a value for the Balmer decrement, assuming 
case B recombination (e.g. Osterbrock 1989; Izotov et al. 1994). 
The detection of several of the Balmer lines in the UVES spec- 
trum, see Table [3] provides a good constraint on the extinction 
value. We assume intrinsic Balmer line ratios at 20000 K, and 
compute flux ratios of the detected Balmer lines. We caution that 
H6 and He are located in the blue end of the UVES spectrum, 
where we see a slight discrepancy in flux calibration between the 
FORS spectrum of March 4 and the UVES spectrum (the UVES 
spectrum having a slightly higher mean continuum flux). We 
find that the Balmer line ratios from the UVES spectrum are all 
(HB/Hy to HB/He) consistent with the theoretical case B recom- 
bination values. The FORS spectra are taken at low resolution 
(R ~ 300), making the extinction derivation based on the Ha/HB 
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Fig. 2. Emission and absorption lines detected in the UVES 
spectrum that were used to derive the host properties. The error 
spectra are plotted with dotted lines. The vertical dashed lines 
denote the mean redshift of the host. The bottom right panel 
shows the [O III] /14364 line from the FORS spectrum. 

ratio from these data less reliable. We test various intrinsic stel- 
lar Balmer absorption strengths, but find no evidence for internal 
reddening from the Balmer decrement in the UVES spectrum, 
with upper limit E(B - V) < 0.03. However this value may be 
influenced by the uncertain flux calibration of the blue end of 
the spectrum, see section [2] The HB and Hy lines are located in 
an area of the spectrum where the FORS and UVES continuum 
fluxes are more in agreement. The ratio of the HB and Hy lines is 
consistent with case B recombination values, E(B - V) < 0.04 (1 
sigma). In Guenther et al. (2006) a reddening of E(B-V) = 0.042 
+ 0.003 is derived from the strengths of the sodium absorption 
lines in the host, and used by Pian et al. (2006) to account for red- 
dening in the host. We choose to use the more common method 
of accounting for extinction by using the Balmer decrement. The 
resulting net extinction (Galactic + host) is close to the value ap- 
plied by Pian et al. (see also Sollerman et al. 2006 for a discus- 
sion) and does not affect further analysis, as the uncertainties in 
the line ratios are dominated by the uncertainties in their fluxes. 
We further note that the high electron temperature seen in this 
source (see section U. 1 ,2I > can make collisional excitation of neu- 
tral hydrogen important, which can mimic reddening and affect 
the Balmer decrement derived reddening. We can not reliably 
evaluate this effect, as we do not measure a significant Balmer 
decrement. 

We use various nebular line flux ratios to evaluate the possi- 
bility of an Active Galactic Nucleus (AGN) contribution to the 
excitation of the nebular lines. Kauffmann et al. (2003) refine the 
popular line ratio diagnostic [O III] ,15008 /HB vs [N II] / Ha, by 
analyzing a large sample of galaxies from the SDSS. They em- 
pirically define the demarcation between starburst galaxies and 
AGN as follows: a galaxy is AGN dominated if 

log([OIII]/15008/H/3) > 0.61/(log([NII]/Ha) -0.05) + 1.3.(1) 

Kauffmann et al. (2003) use extinction corrected fluxes, but note 
that these line ratios are relatively insensitive to extinction ef- 
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fects. We use the flux values measured from the FORS spectra, 
[OIII] /15008/H/? = 4.55 + 0.28 and [Nil] /Ha = 0.06 ± 0.01. 
These flux ratios place the host galaxy comfortably in the lo- 
cus of actively starforming galaxies. Hence we assume that the 
host galaxy nebular emission line excitation is not dominated by 
(non-thermal) AGN emission. 

3.1 .2. Electron density and temperature 

We make the standard two zone model assumption for the H II 
region(s) from which the emission lines originate, consisting of a 
low and intermediate temperature region (e.g. Osterbrock 1974). 
The [Oil] AA3726 and 3729 A lines (Figure© are closely linked 
to collisional excitation and de-excitation at the typical temper- 
atures of star forming regions (~ 10 4 K). The electron density 
(n e ) in the low temperature region can be derived from the ratio 
of the fluxes of these two lines (Osterbrock 1974). The ratio of 
the line fluxes /13726 / /i3729 approaches I in the n e — » limit. 
In the limit n e — > oo the flux ratio approaches 2.86 (collisional 
excitation and de-excitation dominate over radiative transitions, 
forming a Boltzmann population ratio). As shown in Fig. 2, in 
the UVES spectrum the [O II] doublet is clearly resolved, and 
its flux ratio is 0.62 + 0.05. This number is consistent with the 
low density limit, implying that collisional de-excitation is not 
important for the fluxes of the forbidden lines. The resolved 
[Oil] doublet has also been observed in two other GRB host 
galaxy spectra, GRBs 990506 and 000418, which have values of 
[O II] ,13726 / ,13729 of 0.57 ± 0. 14 and 0.75 + 0.11, respectively 
(Bloom et al. 2002). These two detections imply low values for 
n e similar to GRB 0602 1 8 . 

A different doublet used frequently as a density diagnos- 
tic is the [S II] AA6717, 6731 doublet. The line ratio [S II] 
/16717//16731 approaches 1.5 when n e — > 0, and 0.44 above 
n e ~ 10 4 cm 3 . These lines are redshifted out of the UVES 
wavelength range, but are detected in the FORS spectra. We 
find [S II] /16717M6731 = 1.0±0.6, which is therefore not use- 
ful to discriminate between the high and low density regimes. 
Prochaska et al. (2004) find a ratio of 1 . 19 + 0.09 for the host of 
GRB 031203, corresponding to n e ~ 300 cm 4 . We assume the 
relatively low values of n e — 100 cm 4 for our analysis, follow- 
ing e.g. Skillman et al. (1994); Izotov et al. (2006b). 

We estimate the electron temperature (T e ) in the intermediate 
temperature region from the ratio of [O III] nebular and auroral 
fluxes. The auroral [OIII] /14364 is not significantly detected in 
the UVES spectrum - it is located in a noisy region close to the 
gap in between the wavelength ranges, see §2. 

From the FORS spectrum we find the flux ratio [OIII] 
(.14959 + ^5008)/^4364 = 29.1 ± 6.3. The relatively large error 
is mainly due to the rather large uncertainty in the [O III] /14364 
flux value. We use the electron density assumed above, which 
we take to be the same in both the low- and intermediate tem- 
perature regions (see e.g. Osterbrock 1974), and find an electron 
temperature T e (0 2+ ) in the intermediate temperature region of 
2.48 *9'| xlO 4 K. This is high when compared to that of the host 
of GRB 031203 (T e (0 2+ ) ~ 13400 K, Prochaska et al. 2004), as 
shown in Table 1 . Comparably high values of T e (0 2+ ) have been 
observed in the recent discovery of two extremely low metallic- 
ity galaxies by Izotov et al. (2006b). This very high temperature 
and low density suggests a low oxygen abundance for the host 
of GRB 060218, since the main nebular cooling is done through 
oxygen forbidden line emission. 

Especially at low metallicity there can be large differences 
in electron temperature between the low and high temperature 



zones. Due to a lack of detected lines that can be used as tem- 
perature indicators in the low temperature region (the [Oil] 
/L17320, 7331 are redshifted out of the UVES coverage and are 
too faint for the FORS spectrum), we follow the recipe by Izotov 
et al. (2006a), 

T e (0 + ) = -0.577 + T e (o 2+ ) x (2.065 - 0.4987, (o 2+ )) , (2) 

where T e is in units of 10 4 K. We find T e (0 + ) = 1.5+° ' x 10 4 
K, using the approximation for a low metallicity environment 
(12 + log O/H ~ 7.2, see Izotov et al. 2006a). When we as- 
sume an intermediate metallicity (12 + log O/H ~ 7.6), we find 
T e (0 + ) = 1.3+JJ 2 x 10 4 K. We adopt the low metallicity value 
for now, and note that the prescription by Pagel et al. (1992) 
gives a similar value of T e (0 + ) = 1,66+qJ 8 x 10 4 K. The con- 
version above from T e (o 2+ ) to T e (0 + ) is derived through se- 
quences of photoionization models (Stasinska & Izotov 2003). 
Izotov et al. (2006a) show by comparing this conversion with di- 
rect measurements of T e (o 2+ ) and T e (0 + ) that the models agree 
with the measurements. There is a significant scatter of direct 
measurements of T L ,(0 + ) and T e (0 2+ ) with respect to the model 
predictions, which is probably mainly due to the large uncertain- 
ties in the measurements of T e (0 + ) (Izotov et al. 2006a). We do 
not include an additional uncertainty in the following analysis 
to account for this, but note that the small errors reported here 
should not be overinterpreted. We assume the T e (0 + ) to also be 
valid for N II and S II, which have comparable ionisation poten- 
tials. 

3.1.3. Oxygen abundance 

We use the equations from Izotov et al. (2006a) to derive ionic 
abundances from the derived electron densities and tempera- 
tures. To find the oxygen abundance, O/H, we sum the 2+ /H + 
and + /H + ion abundances, assuming the 3+ abundance is 
negligible since no high temperature lines (e.g. He II) were de- 
tected. We take the line fluxes from the UVES spectrum, and find 
2+ /H + = 1.72+0.45xlO~ 5 andO + /H + = 1.78;^ 2 2 xlO- 5 .The 
errors include the uncertainties on the line flux ratios and elec- 
tron temperature. The ratio + / (0 + + 2+ ) > 0.1, confirming 
that the excitation is too low for 3+ to be important, as shown 
by photoionization models by Izotov et al. (2006). The large un- 
certainty in the + /H + abundance is caused by the relatively 
large error in the T e {0 + ). The total oxygen abundance is O/H = 
2+ /H + + + /H + = 3.50 x 10~ 5 , or 12 + log(0/H) = 7.54 
+JJ J 6 , or =* 0.07Z o , assuming log(O/H) + 12 = 8.69 (Allende 
Prieto et al. 2001). This value places the host of GRB 060218 
among the most metal deficient subset of emission line galaxies 
in the local universe (e.g. Izotov et al. 2006a, Lee et al. 2006). It 
is also the lowest metallicity found so far for a GRB host from 
emission line analysis (absorption line metallicities from several 
afterglows show lower line of sight metallicities, down to ~0.01 
Zq for GRB 050730, Starling et al. 2005; Chen et al. 2005). We 
note that in an independent analysis Kewley et al. (2006) found 
a similar value for the metallicity of the host of GRB 060218 of 
12 + log(0/H) ~ 7.6. 

3.2. Relative element abundances 

Accurate emission-line abundances have been derived for a 
small sample of GRB host galaxies. A notable example is the 
spectrum of the host of GRB 031203, for which a solar abun- 
dance pattern was established (Prochaska et al. 2004). We use 
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the host emission lines measured in the UVES and FORS spec- 
tra to gain an insight into the abundance pattern in the host of 
GRB 060218. 

The detection of the forbidden [Ne III] lines allows us to de- 
rive a Ne abundance, using the values for n e and T e (0 2+ ) found 
in §3.1. We use the [Nelll] ,13869/Hy8 flux ratio from the UVES 
spectrum and find Ne 2+ / H + = 2 .9 + 1 .2 x 1 0~ 6 , and Ne 2+ / 2+ = 
0.17 + 0.11, which is consistent with other low metallicity HII 
regions, e.g. in I Zw 18 (where Ne 2+ /0 2+ ~ 0.13, see TableQ}. 
To derive the Ne/H abundance we need an ionisation correction 
function (ICF) for which we use the parametrization by Izotov 
et al. (2006a). We find ICF (Ne 2+ ) = 1.11 andNe/H = 3.3 + 1.3 
xl0~ 6 . 

The [S III] ,16312 line is not detected in the UVES or FORS 
spectra and [S III] AA9532, 9069 are redshifted out of both the 
UVES and FORS coverage. We will therefore only derive a 
value for the ionic S + abundance, and give an upper limit for 
the total sulphur abundance. We use T e (0 + ) and find S + /H + = 
2.7+ 2 {J x 10~ 7 . For the limit on S 2+ /H + we transform T e {0 2+ ) 
to J C (S 2+ ) through the recipe of Izotov et al. (2006a), and use 
the upper limit on the [S III] flux from the UVES spectrum to 
find S 2+ /H + < 2.2 x 10~ 6 . We calculate an ICF(S + + S 2+ ) of 
~1 which allows us to set the not particularly constraining limit 
log(S/0) < -1.1 (Solar value is log(S/0) = -1.50, Lodders 
2003), which is consistent with the observed trend for S to fol- 
low Solar (S/O) ratios independent of (O/H) for low metallicity 
HII regions. 

The [Nil] ,16584 line is redshifted out of the UVES range, 
but is detected in the FORS spectrum. The weaker [N II] ,16548 
line is not significantly detected in the FORS spectrum, with a 
3cr upper limit on the flux of ~ 8x 10~ 17 erg s _1 cm" 2 . We use 
the fixed flux ratio ^6584/^6548 = 2.9 (Osterbrock 1989) and 
T e {0 + ), and find N+/H+ = 1.4+°® x 10~ 6 . 

To calculate N/H from N + /H + we correct for ionisation us- 
ing ICF (N + ) = 1.98, and find N/H = 2.8 + J£ x 10~ 6 , see Figure^ 
The ratio N + /0 + = 0.08 is comparable to, though slightly 
above, the ratio for I Zw 18 of N + /0 + ~ 0.03. We note that 
here we compare the abundances of two elements using two dif- 
ferent spectra (UVES and FORS) taken under different condi- 
tions, and the uncertainty on log(N/0) is likely underestimated. 
Nevertheless the N/O ratio does not significantly deviate from 
the observed trend of low metallicity galaxies to have log(N/0) 
~ -1.5 (e.g. Lopez & Ellison 2003; Izotov et al. 2006). 

Hammer et al. (2006) observed the host of GRB 980425 
/ SN 1998bw with significant spatial resolution (owing to the 
large spatial extent of this host galaxy), and found a high ra- 
tio log(N/0) — -0.6 from a T e abundance analysis at the region 
at which the GRB / SN took place, which corresponds to al- 
most twice the Solar value. This is unexpected at the measured 
metallicity, see Figure[3] Prochaska et al. (2004) find a similarly 
high value of log(N/0) = -0.74 + 0.2 dex in their spectrum 
of the host of GRB 031203, at a metallicity of log(0/H) + 12 
= 8.1. GRB 020903 shows a value more in line with expecta- 
tion from its metallicity (see Fig. 3). In the case of GRB 060218 
the uncertainty on the (N/O) ratio is too high to exclude a de- 
viation from the expected value at the metallicity of the host. A 
possible high N/O value can be explained by a variety of rea- 
sons. If T e has been highly overestimated the metallicity would 
decrease, moving the points in Fig. [3] to the left. Physical rea- 
sons for an enhanced N/O ratio may be e.g. a contribution of 
shock heating to the line emission or a chemical evolution ef- 
fect: Hammer et al. (2006) explain the higher N/O ratio at the 
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Fig. 3. Measurements of log(N/0) and 12 + log(0/H) from a 
sample of galaxies from the SDSS-DR3 (Izotov et al. 2006) are 
shown with open squares, together with points from Lopez & 
Ellison (2003) (small points), for which typical uncertainties are 
indicated in the bottom right corner. The blue dashed lines show 
the tracks where primary and secondary nitrogen dominate the 
N/O ratio. Published GRB host galaxy measurements are over- 
plotted: GRB 060218 (this work), GRB 031203 (Prochaska et 
al. 2004), GRB 980425 (the dashed line connects the value for 
the SN region with the higher log(N/0) and the general host 
value, see §3.2, Hammer et al. 2006) and GRB 020903 (Hammer 
et al. 2006). 



locus of GRB 980425 by a larger N yield of a GRB progenitor 
or SN remnants. 

The hosts of GRBs 980425, 031203, 020903 and 060218 
form a sequence in metallicity, from the metallicities where we 
expect both primary and secondary nitrogen production to play 
an active role (the host of GRB 980425) to where primary N 
is expected to dominate (as in the host of GRB 060218). This 
makes the GRB 060218 host an interesting candidate for deep 
spectroscopy to obtain a more accurate N/O when the SN has 
fully faded. 

The measured nitrogen, oxygen and neon abundances de- 
rived for the host galaxy of GRB 060218 are shown in Table Q] 
We note that these are not spatially resolved. Izotov et al. (1999) 
noted that in the case of I Zw 18 a gradient in electron temper- 
ature can be seen, with the highest temperatures in the regions 
where WR stars are found. These differences in temperature are 
associated with significant differences in (oxygen) abundance 
(with factors up to ~1.4, Izotov et al. 1999). This gradient may 
be due to oxygen enrichment by starforming clusters, and incom- 
plete mixing in the galaxy. Without spatial information we can 
not check abundance gradients in most GRB hosts, and assume 
the oxygen abundances found are representative for the galaxy 
as a whole, including the progenitor locus. However, Hammer 
et al. (2006) and Sollerman et al. (2006) have shown that in the 
case of the host of SN 1998bw strong differences in T e and abun- 
dances are observed as well. 
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-1.574 + 0.06 


Ne 2+ /H + 


2.9 ± 1.2 x 10- 6 




1.91 ±0.16 x 10" 6 


1.83 ± 0.15 x 10-" 


ICF (Ne) 


1.11 




1.20 


1.28 


log (Ne/O) 


-0.8- 


-0.85 ± 0.2 


-0.803 + 0.053 


-0.781 ±0.051 



Table 1. Table of properties of the hosts of GRB 0602 18 (this work) and GRB 031203 (Prochaska et al. 2004) and a comparison to 
the northwest and southeast regions of I Zw 18 (Izotov et al. 1999). 



3.3. Star formation 

The detection of the bright H a /I6563 emission line in the FORS 
spectrum allows us to accurately measure star formation in the 
host: the H a line luminosity only weakly depends on the physi- 
cal conditions of the ionized gas. We use 

SFR H „=7.9xlO- 42 L H „, (3) 

as found by Kennicutt (1998). Moustakas et al. (120061) assess 
the accuracy of this expression through a direct comparison of 
extinction corrected Lh» SFR values and the SFR derived from 
Lir measurements of a sample of IRAS detected galaxies. They 
find that when the H a flux is extinction corrected, the IR and 
Ha SFRs agree without systematic offset with a precision of 
~70%. The extinction corrected flux of Ha can be found from 
the FORS spectra, yielding SFR H „ = 0.065 + 0.005 Moyr 1 . 
Due to slit losses the true H a flux is likely to be higher, and the 
SFRntf can be interpreted as a lower limit. 

Radio and submillimetre observations do not suffer from 
dust extinction. The radio continuum flux of a normal galaxy 
(i.e. non-AGN hosting) is thought to be formed by synchrotron 
emission by accelerated electrons in supernova remnants and by 
free-free emission from HII regions (Condon (1992). It is ex- 
pected that the radio continuum flux is a particularly good tracer 
of the recent SFR, due to the short expected lifetime of the super- 
nova remnants, which is < 10 s yr. We use the method described 
by Vreeswijk et al. (2001) and Berger et al. ( 120031 ) to calculate 
an upper limit to the full star formation rate (ie not influenced 
by any form of dust extinction). We use the deepest 6 cm (4.9 
GHz) Westerbork Synthesis Radio Telescope (WSRT) flux limit, 
i.e. when the initial radio afterglow has faded beyond detection 
limit, and find a 3<x limit of 72 //Jy (formal flux measurement 8 
± 24 //Jy, Kaneko et al. 2006) at a 12h full synthesis on April 
1 2006. This leads to a 3cr SFR upper limit of SFR ladio < 0.15 
Mq yr , which excludes a large amount of obscured star forma- 
tion when compared to SFRh [V . However, we note that despite 
many similarities, GRB hosts studied to date are a diverse popu- 
lation: a handful of hosts show clear indications of much higher 
star formation rates than seen from optical indicators (compa- 
rable to the submillimetre galaxies at several hundred M0 yr -1 ) 
from their submm fluxes (Berger et al. 2003; Barnard et al. 2003) 
and in one case (GRB 980703 at z = 0.97, Berger et al. 2003) its 
radio flux. However, the radio derived SFRs may systematically 
overestimate star formation through a non-negligible contribu- 
tion from AGN activity to the total radio continuum flux. 

One of the properties that sets GRB hosts apart from field 
galaxies is the specific star formation rate (SSFR): the star- 



formation rate per unit mass. This quantity is shown to be high 
in GRB hosts (e.g. Christensen et al. 2004; Courty et al. 2004). 
Sollerman et al. (2006) find that the host of GRB 0602 18 has 
L = 0.008L*, implying a SSFR of ~ 8 M G yr^UV-L.)" 1 - This 
value is comparable to the SSFRs for a sample of GRB hosts 
determined through SED fitting by Christensen et al. (2004). 
Given the multi band SDSS photometry of the host, a low stellar 
mass is expected. After modeling the SED of the host, Savaglio, 
Glazebrook & Le Borgne (in preparation) find M« = 10 7 2±0 3 
Mq. The measured metallicity and stellar mass therefore place 
this galaxy on the mass-metallicity relation found recently in lo- 
cal dwarf galaxies by Lee et al. (2006). 

4. Discrete velocity components in emission and 
absorption 

4. 1 . Emission line profiles 

The resolution of the UVES spectrum is sufficient to search for 
structure in the emission line profiles. Figure |4] shows that the 
brightest emission lines of [O III] AA4959, 5008 significantly de- 
viate from a single Gaussian line profile, and are skewed towards 
the blue. We can rule out an instrumental effect as no such effect 
was seen in the arc or sky line profiles. Due to a lower signal 
to noise, we are not able to verify quantitatively whether other 
emission lines share the same line profile. We fit Gaussian com- 
ponents to the lines, where the number of components and the 
position and width of the components are free parameters. The 
[OIII] AA4959, 5008 lines are fit simultaneously with a fixed 
flux ratio between the /14959 and /i5008 components of 1:3. We 
use both a fit by eye with IRAF splot and a quantitative de- 
blending using vpfitQ. Results agree, and we find a good fit 
(reduced x 2 ~ 0.72) using two Gaussian components, shown in 
Fig.|H where the components are bluewards and redwards of the 
average host galaxy redshift. For the red component we find a 
FWHM of 49 + 5 km s" 1 and redshift 0.033453 + 0.000019; 
and for the blue component a FWHM of 35 + 3 km s~' and red- 
shift 0.033379 ± 0.000005. This gives Az = 7.4 +2.4 x 10~ 5 , 
or ~21.6 km s" 1 velocity separation. These values are similar to 
those found from spatially resolved high resolution spectroscopy 
of emission line regions (i.e. the 30 Doradus nebula, see Melnick 
et al. 1999 for a Har study). Arsenault & Roy (1986) show that 
« 66% of all giant extragalactic HII regions show symmetrical 
Ha lines in their integrated spectra; the remaining profiles show 

1 See http : //www . ast . cam . ac . uk/~rf c/vpf i t . html 
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Fig. 4. The [OIII] /14959 (top) and ,15008 (bottom) emission 
lines. The continuum is normalized. The error spectrum is shown 
by a dotted line. The lines are simultaneously fit using fixed flux 
ratios between ,44959 and /15008 components of 1:3. The two 
seperate Gaussian components are shown with dash-dotted lines, 
and the total emission with a dashed line. 



asymmetries similar to the ones in the host of GRB 060218. 
When asymmetric line profiles are observed (i.e. in high reso- 
lution studies of galactic and extragalactic H II regions) the ex- 
planation of their profiles is usually only possible by using a high 
spatial resolution and correlation of the spectra with images. As 
an example, the integrated Ha line profile of 30 Doradus shows 
a broad and narrow Gaussian component (Melnick et al. 1999), 
whose central wavelengths coincide, while spectra taken at mul- 
tiple positions in the nebula show a large variety of line profiles 
and Gaussian components, which generally can be associated 
with filaments in the nebula. In the case of this host galaxy, we 
observe the integrated profile over 7 kpc. 

There are several possible explanations for the occurence of 
two lines, with the broader one redshifted with respect to the 
blue component. One possibility is that we see two separate star 
forming regions in the host. An alternative explanation is that the 
components are caused by an expanding shell (bubble) around 
the starforming region; or by infalling gas onto the H II region 
(e.g. ionised gas ejected by perhaps SN shocks or stellar winds, 
that falls back). Neither of the last two scenarios seem plausible: 
an expanding shell is not likely to have the measured velocity 
width; and the infalling gas would have to be very highly excited 
to generate the required [O III] luminosity, which makes it diffi- 
cult to have co-exisiting Nal (see §4.2). We tentatively interpret 
the two components as arising from two different starforming re- 
gions in the host. High resolution imaging would be necessary 
to confirm this. There are several GRB hosts where HST images 
resolve the host into multiple starforming regions with similar 
brightnesses (see e.g. Fruchter et al. 2006). We note that this is 
the first identification of resolved emission line components in a 
GRB host galaxy spectrum. 

4.2. Absorption lines in the circumburst medium 

In the UVES spectrum we detect absorption lines at ~ 4070 A, 
~ 4100 A and ~ 6090 A, associated with CaII/U3934, 3969 
and Na I /U589 1,5897 in the foreground gas of the GRB/SN, 
and within the host galaxy. Figure shows at least two discrete 
velocity components, separated in velocity by ~24 km s _1 (sys- 
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Fig. 5. The Nal and Call absorption lines detected in 
GRB 060218. The two velocity components are marked by the 
dashed lines (more precisely, these mark the position of the two 
Na I absorption lines). The smooth line is the result of the best fit 
Voigt profile (reported in Table [2]). The bottom panel shows the 
[OIII] /15008 emission line for comparison. The dotted spectrum 
is the noise. 

tems 1 and 2 in Table |2j see also Guenther et al. 2006). The Ti II 
absorption lines are not detected because they are in a very noisy 
region of the spectrum (A < 3500 A). 

The observed lines have been fitted with Voigt profiles using 
the MIDAS package FITLYMAN (Fontana & Ballester 1995). 
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Table 2. Absorption lines in the UVES spectrum of GRB 060218 



System Line T a i W„ loglV b log (Nal/Call) 
(A) (A) [cm- 2 ] (kms- 1 ) 

1 Na0589l 6088.23 0.033378 0.091 ± 0.008 11.79 ±0.04 6.3 ± 1.0 -0.54 ± 0.08 
Na 1 ,15897 6094.41 0.049 ± 0.007 

Call ,13934 4066.05 0.033363 0.093 ± 0.029 12.33 ±0.07 13.3 ±3.2 
Call ,13969 4102.03 0.064 ± 0.028 

2 Na 1 ,15891 6088.73 0.033462 0.071 ± 0.006 12.22 ± 0.22 1.1 ±0.2 0.01 ±0.23 
Na 1 ,15897 6094.90 0.064 ± 0.006 

Call ,13934 4066.42 0.033458 0.075 ± 0.025 12.21 ±0.08 4.6 ±1.8 

Call ,13969 4102.40 0.054 ± 0.029 

Galaxy Nal ,15891 5891.75 0.0000285 0.33 12.68 ± 0.02 6.4 ±0.2 0.27 ±0.06 

Na 1 ,15897 5897.73 0.25 

Call ,13934 3934.93 0.0000381 0.13 12.41 ±0.06 (5) 

Call ,13969 3969.74 0.08 
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Fig. 6. Na I and Ca II column densities for the two systems in the 
host of GRB 060218 (filled dots) and lines of sight in the Milky 
Way (open squares; typical error for these cases is < 0.1 dex; 
Hunter et al. 2006). 



Figure [5] shows the fitting results, and column densities and 
Doppler parameters are listed in Table [2] The lines are barely 
resolved, i.e. the line widths are at the level of the UVES spec- 
tral resolution FWHM ^6.5 km s~' . 

The measured FWHM corresponds to an instrumental PSF 
of 3.9 km s" 1 , if expressed in terms of Doppler parameter, which 
can be translated into an upper limit on the gas temperature (de- 
rived from the lightest element Na, assuming thermal broaden- 
ing) of T < 2 X 10 4 K. Narrow metal lines are often detected in 
GRB afterglows when high resolution spectra are acquired (see 
for instance Chen, Prochaska, & Bloom 2006). 

There are indications that Na I and Ca II are not tracing each 
other. The positions of the two Ca II components are shifted by a 
few km s~' with respect to Na I, which could be consistent with 
uncertainties on the fitted parameters (more severe for the Ca II 
doublet). However, the line broadenings and Ca II/Na I ratios are 
different in the two components. We note that in the Galactic 
ISM Call and Nal are found in regions with different physi- 
cal conditions. The former is found to be generally broader than 



the latter, indicating that it traces warmer, more turbulent, and/or 
larger gas clouds (Welty et al. 1996). 

Remarkably, the positions of the two absorption systems are 
consistent with the redshifts of the two emission-line compo- 
nents in the HII regions derived independently (see §3.5 and 
the lower panel of Figure[5]l. The difference in redshift between 
emission and absorption is < 1 km s" 1 and ~ 2 km s~' for the 
blueshifted and redshifted systems, respectively. The relative ve- 
locity for the two emission components is 21 km s _1 , close to 
the 24 km s~' measured from the absorption systems. However, 
the broader component in absorption is at the lowest redshift, 
whereas the opposite is true for the emission. 

The Ca II column density was measured in another two GRB 
afterglows (Savaglio & Fall 2004; Savaglio 2006) with a total 
column density of nearly 10 14 cirT 2 in each of them. Na I absorp- 
tion in GRBs is reported here for the first time, basically due to a 
lack of suitable data in past GRB observations (Na I is redshifted 
into the NIR for z > 0.7). The Nal and Call abundances have 
been studied in the Galaxy and LMC (Hunter et al. 2006; Vladilo 
et al. 1993). Beyond the Local Group, Nal and Call have been 
detected in 2 damped Lyman-o' systems (DLAs) in QSO spectra, 
at z = 1.062 and 1.181 (Petitjean et al. 2000; Kondo et al. 2006). 
Other QSO absorption line studies report upper limits for Nal 
(Boksenberg et al. 1978; 1980). GRB 060218 is the third source 
outside the Local Group where Nal is measured^ 

It is rather complicated to interpret the detection of Na I and 
Call in GRB 060218 in terms of relative abundances. The ion- 
ization potentials of the two ions are quite different (5.1 eV and 
1 1 .9 eV for Na I and Ca II, respectively). Na II is likely the dom- 
inant ion in a neutral-gas environment (the H I ionization poten- 
tial is 13.6 eV), whilst Calll dominates the calcium species. To 
derive the total abundance of calcium, a significant ionization 
correction can therefore be necessary. Moreover, Ca is an a el- 
ement, while Na is not. Hence in low metallicity systems, like 
GRB 060218, an a-element enhancement is expected. Sodium 
was found to trace iron in stars with metallicities close to Solar, 
but it can have lower abundances for lower metallicities (Timmes 
et al. 1995). Our expectations are further complicated by the 
rather different refractory properties of the two elements: Na is 
little depleted on to dust grains, whereas Ca can be 99% depleted 
(Welty et al. 1994; Savage & Sembach 1996). This problem may 
be somewhat mitigated by the likely negligible dust depletion in 



2 More Nal lines are detected in low z SDSS QSO spectra, but no 
column density measurements are reported for these systems (V. Wild, 
private communication.) 
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the gas, as suggested by the small dust extinction derived in the 
H II regions of GRB 060218 (see §3.1). 

Nevertheless, we compare Nal and Call in GRB 060218 
with what is typically observed in the ISM of the Milky Way 
(Figure |6}. The two systems in the host of GRB 060218 lie in 
the bottom left corner of the distribution. Is the observed Na and 
Ca behaving like the neutral gas of the Milky Way? If the ab- 
sorption lines are arising in a neutral region of the ISM, and if 
we consider the empirical relation that links H I to Na I (derived 
by Hunter et al. 2006) we would expect an H I column density 
along the GRB sight line of the order of log/Vei = 20.6 cm -2 . 
However, if the metallicity in the neutral gas is similar to the 
H II regions probed by the emission lines, a much lower total H I 
column density is more likely: logA^i ~ 19.4, suggesting that 
perhaps the Na I-H I relation found for the Milky Way may not 
be applicable for the gas around GRB 060218 or that the metal- 
licity in the galaxy is not uniform. An H I column density of the 
order of 10 19 4 cm~ 2 is relatively large, considering that the stel- 
lar mass of this GRB host is more than 1000 times smaller than 
that of the Milky Way. However, a large reservoir of gas is ex- 
pected given the high SFR per unit stellar mass estimated for the 
host (see section [331 ). 

5. Secondary metallicity calibrators 

The association of GRB 060218 with a low mass, low metal- 
licity, high excitation host galaxy follows the trend seen from 
the other nearby GRBs (e.g. Sollerman et al. 2005). In fact, the 
very low metallicity of the host significantly extends the known 
metallicities of GRB hosts through emission line spectroscopy. 
If the redshift of the galaxy had been substantially higher, the 
[OIII] ,14364 line would have gone undetected, making a di- 
rect determination of the abundance through T e impossible, and 
the metallicity inferred from secondary calibrators would have 
placed the host at significantly higher metallicity (e.g. Mirabal 
et al. 2006; Modjaz et al. 2006). The sample of GRB host galax- 
ies with a measure of abundances through T e is limited to just 
four nearby sources due to the faintness of the [OIII] ,14364 
line, where GRB 020903 has the highest redshift with z = 0.25. 
For all other galaxies we are forced to use secondary, empiri- 
cal methods to calculate metallicity. The most common is the 
R23 calibrator (see e.g. Kobulnicky & Kewley 2004), which 
uses the bright [OIII], H/3 and [Oil] lines. This method pro- 
duces a degenerate metallicity solution, which can be broken 
through other emission lines (e.g. [Nil] and Ha), but in many 
cases these other lines are not available (due to e.g. redshift or 
insufficient S/N). The R23 method is calibrated through photo- 
ionization models, which have limitations at the low and high 
metallicity ends. It has been shown that R23 metallicities have 
an offset with respect to T e metallicities for metallicities close 
to Solar (e.g. Bresolin et al. 2004). Modjaz et al. (2006) use the 
R23 method to find 12 + log(0/H) = 8.0 + 0.1 for the host of 
GRB 060218. An alternative method that is calibrated on a sam- 
ple of H II regions with T e determined metallicities, is the ra- 
tio of the nebular [O III] and [N II] lines (e.g. Pettini & Pagel 
2004). The N2 index (N2 = log([N II],16583/Ha)) has been 
proposed particularly for low metallicity galaxies. However, for 
GRB 060218, the N2 calibration overestimates the metallicity 
of the host by more than a factor of 3, yielding 12 + log(0/H) ~ 
8.2, but the scatter in this relation is large (90% of the fitted data 
in Pettini & Pagel 2004 falls within a range log(0/H) = +0.4). 
The 03N2 = log (([O III]^5008/H/3)/([N 1 1] ,165 8 3 /Ha)) ratio 
has significantly less scatter but above 03N2 > 1.9 this indi- 
cator breaks down and cannot be reliably used (Pettini & Pagel 



2004). Our host has 03N2 = 1.88 + 0.11, which may explain 
the overestimate of the metallicity through this indicator of 12 + 
log(0/H)^8.1. 

Concluding, we can state that an analysis using a variety of 
different secondary methods (e.g. R23, 03N2) would not have 
found the true metallicity, but may still have identified this host 
as a very low metallicity candidate. 

6. Massive stars and progenitors 

Long GRB progenitors are likely massive Wolf-Rayet (WR) 
stars. We can gain a greater understanding of the evolutionary 
paths of such massive stars towards GRBs by detecting WR pop- 
ulations within GRB host galaxies. A valuable diagnostic on the 
massive star population is the He II /14686 line, which appears 
as a broad line and/or as a nebular line in some HII galaxies 
(Schaerer et al. 1999), and is a direct sign of (unusually) high 
excitation levels caused by the presence of WR stars (especially 
WC- and WO-type WR stars). Together with, amongst others, 
N III ,14640, C IV ,14658, [Fe III] ,14658 and [ArlV] ,1471 1 lines, 
this line forms the so-called blue WR bump (from ~4650 - 4700 
A) in low resolution spectra, which is often accompanied by a 
C IV ,15808 line (the red bump). In a recent deep spectroscopic 
search in nearby GRB host galaxies, Hammer et al. (2006) have 
detected the He II line and accompanying WR bump in the spec- 
tra of the hosts of GRB 980425, 020903 and 031203. From the 
relative intensities of the H /3 and He II /14686 lines or WR bump 
(measured flux ratios H (5/ He 11 ,14686 generally range from 
~0. 01-0. 02) one can estimate the ratio of WR to O stars. 

In Crowther & Hadfield (2005) the effect of metallicity on 
the derived WR / (WR + O) ratio is investigated, by comparing 
the fluxes of SMC, LMC and Galactic WR stars with atmosphere 
models. They find a lower WR line luminosity at decreasing 
metallicity. In the hosts of GRBs 980425 and 020903, Hammer 
et al. (2006) find values of WR/O ~ 0.05 and ~ 0.14 - 0.2, 
respectively. We note that no metallicity correction has been ap- 
plied to these values, which would increase the number of WR 
stars as the studied hosts have sub-Galactic metallicity. 

These high WR / O ratios are of particular interest as their 
metallicities are accurately known from T e analyses: 0.5 Zq 
and 0.19 Z for GRB 980425 and GRB 020903, respectively 
(Hammer et al. 2006). At these low metallicities the high mea- 
sured WR / O ratios may be regarded as evidence for very recent 
star bursts, but it is difficult to determine the ages of the dom- 
inant stellar populations. Optical and near-infrared SED fitting 
has shown that the dominant stellar populations in GRB hosts 
are young at typically ~0.1 Gyr (Christensen et al. 2004). In the 
case of the host of GRB 980425 the entire galaxy is well fitted 
with a continuous star formation history (Sollerman et al. 2005), 
although Hammer et al. (2006) find a young population from 
the equivalent widths of the emission lines. No strong evidence 
for a very young starburst is apparent in the host spectrum of 
GRB 020903, but the merger morphology of the host and the 
emission line EWs may suggest that some recent star formation 
is present. However, the WR/O ratio of 0.14 - 0.2 in the host of 
GRB 020903 is particularly remarkable, as such abundant pro- 
duction of WR stars at the low metallicity of 0.19 Zq can only 
be explained by invoking instant star bursts with peculiar initial 
mass function within the standard star burst model by Schaerer 
& Vacca (1998) (e.g. Fernandes et al. 2004). 

Recent stellar evolution models indicate that the effects of 
rotation may be, in part, responsible for observed high WR / O 
ratios in galaxies. According to Meynet & Maeder (2005a), in- 
cluding the effects of rotation significantly enhances the mass 
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loss rates of massive stars during the giant phase compared to the 
non-rotating case, as the shear instability due to the strong de- 
gree of differential rotation between the core and the envelope in- 
duces fast chemical mixing. Their models predict a WR/ O ratio 
of about 0.02 at Z = 0.004 for a constant star formation, and in 
principle WR/O ~ 0.2 might be achieved at the given metallic- 
ity for instant star bursts even with a standard initial mass func- 
tion. However, such rotating models predict GRB / SN Ibc ratios 
that are too high. Their prediction of spin rates of young neutron 
stars is also inconsistent with observations (Hirschi, Meynet & 
Maeder 2005; see also Heger, Langer & Woosley 2000). 

More recent models that include magnetic torques for the 
angular momentum transport in the star suggest another way to 
produce WR stars at low metallicity. Although in magnetic mod- 
els the chemical mixing induced by the shear instability during 
the giant phase is negligible as the magnetic torque tends to keep 
the star rigidly rotating, the mixing by meridional circulation 
can be very fast even on the main sequence (Maeder & Meynet 
2005b). This may even cause the whole star to become chemi- 
cally homogeneous if the initial rotational velocity is sufficiently 
high, especially at sub-solar metallicity (Yoon & Langer 2005; 
Woosley & Heger 2006). 

Formation of WR stars can thus be induced not only by 
mass loss but also by chemical mixing, and Woosley & Heger 
(2006) found that some WR stars formed through such chemi- 
cally homogeneous evolution can actually retain enough angu- 
lar momentum necessary for GRB production. The predicted 
GRB / SN ratio through such evolutionary channels also turns 
out to be consistent with observations, when the observationally 
derived initial rotational velocity distribution of massive stars 
by Mokiem et al. (2006) is adopted (Yoon, Langer & Norman 
2006). Interestingly, as stars may be kept rotating rapidly at low 
metallicity due to lowered mass loss rates (e.g. Vink et al. 2005), 
this new type of evolution by fast rotational mixing can lead to 
high WR / O ratios even at very low metallicity (Yoon, Langer & 
Norman 2006). It also predicts rather large delay times for WR 
production from the star formation (i.e., from several 10 6 yrs to 
a few ~ 10 7 yrs; Yoon, Langer & Norman 2006) compared to 
those from mass-loss induced WR formation ( < a few 10 6 yrs). 

In this regard, an estimate of the WR / O ratio in the host 
galaxy of GRB 060218 might provide a valuable test case for 
the new GRB progenitor scenario by Yoon & Langer (2005) and 
Woosley & Heger (2006). That is, if the WR / O ratio in this host 
turns out to be high despite its very low metallicity, it may sup- 
port the chemically homogeneous evolution scenario for GRB 
progenitors and extend the metallicity - WR/O star ratio rela- 
tion for GRB hosts down to lower metallicity. This relation is 
well known for the Milky Way, LMC and SMC, and a devia- 
tion of GRB hosts away from that trend gives strong input for 
progenitor modelling. 

The high resolution UVES spectrum of GRB 060218 should 
be able to resolve the components comprising the WR bump, 
and resolve a He II line into a broad and a nebular component. 
However, no nebular He II line or WR bump has been signif- 
icantly detected. We measure the flux upper limit on the WR 
bump in the UVES spectrum by summing the flux in the WR 
bump wavelength region (restframe 4650 - 4686A). Following 
the method of Schaerer & Vacca (1998) we set an upper limit of 
WR / (WR + O) < 0.4, which is not a constraining limit, owing to 
the fact that the SN outshines the possible He II line. To reliably 
detect the WR bump in this host we need to detect the host con- 
tinuum with reasonable S/N, which means the supernova has to 
fade below this level before more sensitive searches are feasible. 



7. Conclusions 

We present a VLT UVES high resolution spectrum of SN 2006aj, 
associated with the nearby GRB 060218 at heliocentric redshift 
z = 0.03342(2). We use the emission lines of the UVES spec- 
trum as well as the line measurements from our FORS spectro- 
scopic campaign to derive properties of the host. We find that 
the electron density is low, and the electron temperature is high, 
T e (0 2+ ) = 2.48+^ x 10 4 K, as shown in Table 1. We find a low 
host metallicity of 12 + log(0/H) = 7.54+q] 7 ,, placing it among 
the most metal deficient subset of emission-line galaxies. It is 
also the lowest metallicity found so far for a GRB host from an 
emission line analysis. The metallicity we find lies considerably 
below the values derived using secondary calibrators, e.g. the 
metallicity 12 + log(0/H) = 8.0 as derived from the R23 calibra- 
tor (e.g. Modjaz et al. 2006). The mass of the galaxy is low, and 
matches what is expected from the mass-metallicity relation for 
dwarf galaxies. We measure a relatively high value for log(N/0) 
with respect to the metallicity, which is also seen in a few other 
GRB hosts. As our uncertainty on log(N/0) is relatively high, 
deeper spectroscopy is needed to confirm this overabundance. 

The bright emission lines show strong evidence for asymme- 
try, and a single Gaussian provides a poor fit to the profiles of the 
bright [OIII] emission lines. A two Gaussian model provides a 
satisfactory fit with the two components separated by ~22 km 
s _1 . We find the same two velocity components in absorption 
through the Ca II and Na I absorption lines in the host. We tenta- 
tively interpret these two velocity components to be due to two 
star forming regions in the host galaxy. However, to unravel their 
true identity, high spatial resolution imaging is needed. 

The dust content of the galaxy is low, based on the Balmer 
line decrement. This is also evident from the low limit on the 
obscured star formation rate we set through a 3<x upper limit on 
the flux at 6 cm of SFR ra d; < 0.15 M0yr _1 , compared to the 
optical star-formation rate SFR Hff = 0.065 ± 0.005 M Q yr _1 . 

This host galaxy is an interesting target for future spec- 
troscopy targeted at the WR bump, as the low metallicity of the 
host will significantly extend the present sample of GRB hosts 
with known WR star content and metallicity. We show that a 
measure of these two quantities for a sample of GRB hosts may 
provide further insight into the nature of GRB progenitors. 

The absolute magnitude of the host (Mb = -15.9, 
e.g. Sollerman et al. 2006) is such that this galaxy would not 
have been detected in any survey at a redshift of z ~ 1, let 
alone at the mean Swift GRB redshift of z ~ 2.8 (Jakobsson 
et al. 2006), which makes this host an important object to study 
in the context of larger redshift GRB hosts. 
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